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Scalar  relativistic  self  consistent  linearized  muffin-tin  orbital  (LMTO) 
calculations  for  bcc  and  fee  Fe  have  been  performed  with  several  dif¬ 
ferent  local  approximations  to  the  exchange  and  correlation  energy 
density  and  potential.  Overall,  in  contrast  to  the  conclusions  of  previ¬ 
ous  studies,  wo  find -eh at  the  local  spin  density  approximation  (LSD A) 
to  exchange  and  correlation  can  provide  an  adequate  description  of 
bulk  Fe.  provided  a  proper  parameterization  of  the  correlation  energy 
density  and  potential  of  the  homogeneous  electron  gas  over  both  spin 
and  density  is  used.  Lattice  constants,  found  from  the  position  of 
the  minimum  of  the  total  energy  as  a  function  of  Wigner-Seitz  radius, 
agree  to  within45&4ior  .s,ptdLMTOs  onlyj  and  withfft  1-2%,  (for  s,p.d,f  . 
LMTOs^of  the  experimental  lattice  constants  For  all  forms  used  for' 
the  local  correlation.  The  best  agreement,  however,  was  obtained  us¬ 
ing  a  local  correlation  potential  derived  from  the  Vosko-Wilk-Nusair 
form  for  the  spin  dependence  of  the  correlation  energy  density.  The 
calculation  performed  with  this  correlation  potential  was  also  the  only 
calculation  to  correctly  predict  a  bcc  ferromagnetic  ground  state.  All 
the  calculations  gave  satisfactory  values  for  the  magnetic  moment, 
though  the  experimental  bulk  modulus  was  overestimated  by  more 
than  30%.  a  feature  in  agreement  with  most  previously  published 
calculations.  NThe  agreement  in  bulk  modulus  is  comparable  to  that 
obtained  in  paramagnetic  calculations;  therefore  we  attribute  this  er¬ 
ror  not  to  a,  failure  in  the  LSDA.  but  rather,  this  error  is  most  likely 
a  result  toahe  atomic  sphere  approximation  (ASA). 
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Scalar  relativistic  self  consistent  linearized  muffin-tin  orbital  (LMTO)  calcula¬ 
tions  for  bcc  and  fee  Fe  have  been  performed  with  several  different  local  approxi¬ 
mations  to  the  exchange  and  correlation  energy  density  and  potential.  Overall,  in 
contrast  to  the  conclusions  of  previous  studies,  we  find  that  the  local  spin  density 
approximation  (LSDA)  to  exchange  and  correlation  can  provide  an  adequate  de¬ 
scription  of  bulk  Fe,  provided  a  proper  parameterization  of  the  correlation  energy 
density  and  potential  of  the  homogeneous  electron  gas  over  both  spin  and  density 
is  used.  Lattice  constants,  found  from  the  position  of  the  minimum  of  the  total 
energy  as  a  function  of  Wigner- Seitz  radius,  agree  to  within  1%  (for  s,p.d  LMTOs 
only)  and  within  1-2%  (for  s,p,d,f  LMTOs)  of  the  experimental  lattice  constants  for 
all  forms  used  for  the  local  correlation.  The  best  agreement,  however,  was  obtained 
using  a  local  correlation  potential  derived  from  the  Vosko-Wilk-Nusair  form  for  the 
spin  dependence  of  the  correlation  energy  density.  The  calculation  performed  with 
this  correlation  potential  was  also  the  only  calculation  to  correctly  predict  a  bcc 


ferromagnetic  ground  state.  All  the  calculations  gave  satisfactory  values  for  the 
magnetic  moment,  though  the  experimental  bulk  modulus  was  overestimated  by 
more  than  30%,  a  feature  in  agreement  with  most  previously  published  calculations. 
The  agreement  in  bulk  modulus  is  comparable  to  that  obtained  in  paramagnetic 
calculations;  therefore  we  attribute  this  error  not  to  a  failure  in  the  LSDA,  but 
rather,  this  error  is  most  likely  a  result  of  the  atomic  sphere  approximation  (ASA). 

PACS  Indices:  75.50.Bb,  71.10.-f-x 
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INTRODUCTION 


Considerable  success  has  been  achieved  by  electronic  structure  calculations  in  the  pre¬ 
diction  of  the  ground  state  properties  of  non-magnetic  bulk  crystals.  Lattice  constants  can 
usually  be  obtained  to  within  1  —  2%  or  so  of  the  experimental  value}  elastic  constants 
such  as  bulk  moduli  to  within  0-30%  }'2  Consistent  results  have  been  obtained  using  a 
variety  of  electronic  structure  techniques,  including  the  linearized  augmented  plane  wave 
approach  (LAPW)  (or  its  generalization  to  full  potential,  the  FLAPW  method'*).  Korringa- 
Kohn-Rostoker  (KKR)1,  linearized  muffin- tin  orbital  (LMTO)4,  full  potential  LMTO5  and 
pseudopotential6  approaches  to  name  a  few.  All  of  these  methods  adopt  the  local  den¬ 
sity  approximation  (LDA)  to  the  exchange- correlation  energy  and  potential.  There  have 
been  many  forms  suggested  for  the  energy  and  potential  functionals.  The  commonly  used 
classes  are:  (1)  the  form  suggested  by  Slater7  (Xat)  and  the  exchange  only  Kohn-Sham 
(KS)8  form  where  the  functional  dependence  on  the  density  is  of  the  simple  ps  form,  (2) 
the  fits  to  RPA  based  calculations  by  Hedin  and  Lundqvist  (HL)?  Gunnarsson,  Lundqvist, 
and  Wilkins  (GLW),10  and  von  Barth-Hedin  (VBH),11  and  (3)  the  parameterization  of  the 
Monte  Carlo  data  of  Ceperley  and  Alder12  by  Perdew  and  Zunger  (CA)13  The  latter  is 
expected  to  be  the  best  test  of  the  LSDA,  since  it  provides  the  most  accurate  description 
of  the  homogeneous  electron  gas. 

Historically  the  LSDA  has  not  had  the  same  success  in  magnetic  systems  as  the  LDA 
has  had  in  non-magnetic  systems.  The  3d  transits  .  metal  Fe  is  a  good  example  of  a 
magnetic  system  which  has  been  extensively  studied,  yet  one  which  has  been  to  date 
rather  poorly  described  by  LSDA-based  band  structure  calculations.  The  structural  and 
magnetic  properties  of  Fe  have  been  studied  using  the  FLAPW  method  by  Wang  et  a J.14 
Hathaway  et  a I.15  and  Jansen  et  ai.16  In  these  calculations,  the  authors  find  that  the 
ground  state  crystal  structure  of  Fe  is  non-magnetic  fee,  that  the  bulk  modulus  is  too 
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large  and  that  the  lattice  constant  is  significantly  too  small.  The  original  calculations 
of  Kubler17  also  produced  an  incorrect  ground  state.  More  recent  LMTO  calculations  by 
Bagno  et  al. 19  produced  similar  results  to  the  FLAPW  studies  though  a  good  bulk  modulus 
was  obtained,  unlike  any  other  LSDA  calculation.  The  authors  also  considered  the  effects 
of  various  gradient  terms  (non-local  corrections)  where  they  found  that  they  could  indeed 
stabilize  the  ferromagnetic  bcc  phase  and  also  expand  the  lattice  constant.  It  is  the  purpose 
of  this  Rapid  Communication  to  consider  the  sensitivity  of  the  aforementioned  ground  state 
properties  to  changes  in  the  form  of  spin  dependence  of  the  correlation  effects,  perhaps 
answering  the  questions  of  why  previous  LSDA  calculations  have  given  a  poor  ground  state 
description  and  whether  it  is  necessary  to  resort  to  non-local  corrections. 

All  the  LSDA  calculations  are  based  upon  fits  of  the  exchange-correlation  energy  den¬ 
sity  over  the  electron  gas  density  (or  equivalently,  rs)  and  spin  polarization,  (.  The  latter 
is  usually  approximated  at  intermediate  (  by  an  interpolation  between  the  extreme  para¬ 
magnetic  <(,  =  0)  and  ferromagnetic  (£  =  1)  limits  using  a  formula  suggested  by  von  Barth 
and  Hedin.11  This  scaling,  while  correct  for  the  exchange  part  of  the  potential  and  energy 
in  a  homogeneous  system,  is  not  correct  for  the  correlation  contribution.  Vosko  et  ai.1^ 
(VWN)  suggested  an  improved  spin  interpolation  formula  for  the  correlation  energy,  which 
successfully  fits  the  CA  data  over  the  range  of  (  to  within  1  mRy,  and  it  gives  a  significant 
improvement  over  the  VBH  interpolation.  This  interpolation  for  the  correlation  energy  is 
of  complex  form,  and  the  corresponding  functional  forms  for  the  potentials  are  therefore 
also  complex,  since  the  potentials  arc-  derived  from  the  functional  derivatives  of  the  energy. 
The  equations  for  the  correlation  energy  and  potentials  sis  a  function  of  spin  polarization 
are  as  follows:20 


«c(r„C)  =  <f(r«)  +  A«c(r„C) 
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(1) 


(2) 


where  A ec(rs,()  is  expressed  in  terms  of  the  spin  stiffness  ae 

Aec(r„C  )  =  ae(r,)HN  + /?(r,)C4] 


and 


0  -  4 + A. + m.  [-^  (1 + «,,<<) + ««(.)  (-j^c4) 


3  drs 

[4«rs)(3/«)  +  (l  +  3(r,X4)  /(C)]  (1  =F  0 


where 


0(r$)  = 


/"(0)Aec(r„l) 


Mra) 


-  1, 


^(rs)  =  «f(r«)  ~ 


rs 


def(rs) 


3  dr« 


and 


/(C)  =  [(1  +  C)3  +  (1  -  <)*  -  2]/(25  -  2). 


(3) 

(4) 

(5) 

(6) 


Vosko  et  al.19  obtained  the  following  expression  for  ac 
ac(r8)  =  Aq  ^ln  Js-  ;  +  ^ctan-1  ^ 


_  _ _ bxp  r,  (v^-sp)2 

*(rs)  ^  C?""  2^+6  *Y(iq)  X(rs) 


,  2(6  +  2i0).__-i 

H - - - tan 


Q 


2,/rJ  +  b 


0 
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where 


-Y(rs)  =  r,  +  6>/r7  +  c;  Q  =  y/ Ac  -  b2 


(8) 


and  where  appropriate  parameters  are  given  in  Table  1.  The  derivative  of  ac  is  given  by 


dac(rs)  . 

r,  3 - =  Aa 

dr  g 


1  4-  b\rs? 


1  2 
,1  +  b\rsi  +  &2rs  +  &3rs2 


(9) 


where  &i  =  62  =  and  63  =  - 


cz  0 ' 


CXO  ’  *  CZO 

The  derivative  of  0  is  given  by 

d0(r,)  /"(  0)  d(ef  (*"«)  “  «e*  (r«))  /"(0)  (ef  (r,)  -  dac(rs) 

"  '  a2(rf)  "  dr. 


dr,  ac(r,) 


dr. 
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In  the  above,  P  denotes  the  paramagnetic  limit  and  F  the  ferromagnetic  limit.  In  these 
limits,  we  have  used  the  Perdew-Zunger  13  parameterization: 


/  l+fijh+fr.'  r*~l  (11) 

l  .4;lnrs  +  B{  +  C,r«lnrs  +  Z?,rs,  rs  <  1 


(12) 


and  ,  .  , 

f  5i(iii£nAl±Mri}  _  >  i 

n'c  =  {  l+Pl^+ftr.  ’  4  - 

[  i4jlnrs  +  (B{  -  3A,)  +  |Ctrslnrs  +  j(2 D,  -  Ct)rs,  rs  <  1 

The  VWN  parameterization  for  and  ,  which  is  distinct  from  the  VWN  spin  interpola¬ 
tion,  could  have  also  been  used.  Although  this  makes  comparison  with  other  VWN  results 
difficult,  the  difference  between  the  two  is  less  than  1%  19  over  all  physical  ranges  of  rs. 
The  results  here  were  intended  to  emphasize  the  sensitivity  of  ground  state  properties  to 
the  form  of  spin  interpolation. 


The  values  of  the  various  parameters  are  given  in  table  2,  and  all  energies  are  in 
Rydbergs. 

The  failure  in  previous  calculations  using  the  LSDA  then  could  arise  from  either  the 
breakdown  of  the  LSDA  itself  or  from  the  use  of  the  VBH  spin  interpolation  formula  in 
obtaining  the  spin  dependent  correlation  potential.  We  tested  a  number  of  parameterized 
LSDA  functionals  to  see  whether  the  relatively  poor  agreement  with  experiment  obtained 
in  previous  calculations  for  magnetic  systems  compared  to  non-magnetic  systems  can  be 
remedied  by  using  another  form  for  the  spin  dependent  correlation  potential. 
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CALCULATIONAL  DETAILS 


The  calculations  presented  in  this  paper  were  performed  with  the  scalar  relativistic 
LMTO  method,  which  used  a  fully  relativistic  frozen  core.  We  have  tested  several  different 
local  density  functionals:  namely  the  VBH,  CA,  KS,  and  GLW  forms  with  the  VBH 
spin  interpolation  and  the  CA  with  the  VWN  interpolation.  We  found  that  for  all  the 
functionals  tested  (CA,  VBH,  and  GLW)  that  the  bcc  ferromagnetic  phase  had  a  lower 
total  energy  than  the  corresponding  non-magnetic  bcc  phase,  though  we  did  not  test  the 
stability  of  antiferromagnetic  fee  phase  relative  to  the  non-magnetic  phase.  Care  was  taken 
both  to  converge  the  total  energies  to  6  decimal  places  and  to  sample  the  region  around  the 
minimum  densely  so  that  the  lattice  constant,  magnetic  moment,  and  particularly  the  bulk 
modulus  could  be  accurately  determined.  The  number  of  k  points  used  in  the  Brillouin  zone 
integration  was  systematically  increased  until  convergence  was  achieved.  The  tetrahedron 
method  was  employed  for  this  integration  with  convergence  typically  achieved  by  500  k 
points.  LMTO  calculations  with  s,p,  and  d  as  well  as  s,p,d,  and  f  LMTOs  have  been 
examined.  For  non-magnetic  fee  Fe  and  ferromagnetic  bcc  Fe,  we  summarize  the  results  of 
calculations  based  on  1=3  with  about  500  k  points  for  the  CA  and  VBH  functionals  only. 
The  same  trends  were  also  observed  in  the  1=2  calculations;  the  addition  of  the  f  LMTO 
contracted  the  lattice  and  reduced  the  moment  slightly.  The  details  and  approximations 
behind  the  LMTO  approach  have  been  discussed  in  detail  in  other  works,  and  we  refer  the 
reader  to  reference  4  for  an  excellent  discussion  of  this. 
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A 


RESULTS 


Figure  1  shows  results  of  the  energy  versus  Wigner-Seitz  radius  for  non-magnetic  fee 
Fe  obtained  using  the  VBH  and  CA  density  functionals  respectively  for  525  k  points.  The 
total  energy  plotted  excludes  terms  from  the  frozen  core.  Since  we  are  always  interested 
in  energy  differences,  this  constant  will  not  be  important.  Figure  2  shows  calculations 
with  506  k  points  for  ferromagnetic  bcc  Fe  with  the  VBH  functional  and  the  CA  func¬ 
tional  using  both  the  VBH  and  VWN  spin  interpolation  formulae.  The  results  for  lattice 
constant,  magnetic  moment,  bulk  modulus,  and  stable  ground  state  crystal  structure  are 
given  in  table  3.  We  see  that  lattice  constants  are  within  2%,  magnetic  moments  within 
5%,  and  bulk  moduli  are  too  high  by  30-40%.  The  errors  in  bulk  moduli  are  comparable 
to  those  obtained  in  paramagnetic  LMTO  calculations  for  the  3d  transition  metals.  The 
bulk  modulus  is  extremely  sensitive  to  the  curvature  of  the  total  energy  curve,  and  errors 
associated  with  the  LMTO  approximation  aside  from  LSDA  could  easily  account  for  this 
discrepancy.  Work  to  elucidate  the  source  of  error  is  at  present  underway  using  the  full 
potential  LMTO  method?1  The  magnetic  moment  is  slightly  underestimated.  This  is  a 
result  of  the  Wigner-Seitz  sphere  approximation  which  will  tend  to  overstimate  the  contri¬ 
bution  from  the  antiferromagnetically  polarized  interstitial  region.  The  significant  success 
of  the  new  parameterized  correlation  potential  functional  lies  in  the  correct  prediction  of  a 
ferromagnetic  bcc  ground  state  -  a  feature  not  observed  in  those  functionals  which  use  the 
VBH  form  of  spin  interpolation.  The  improved  VWN  spin  interpolation  in  the  correlation 
potential  is  consequently  important  for  an  accurate  description  of  the  magnetism. 
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CONCLUSIONS 


Our  results  show  that  Fe  has  a  sufficient  sensitivity  to  correlation  that  the  low  lying 
states  can  reorder  with  the  application  of  different  spin  interpolation  forms.  We  used  the 
spin  interpolation  formula  of  Vosko  et  ai.1^  for  the  correlation  energy  density  to  provide  ac¬ 
curate  parameterized  correlation  potentials  within  the  LSDA.  The  VWN  spin  interpolated 
correlation  potential  significantly  improves  the  LSD  description  of  ferromagnetic  bcc  Fe. 
While  the  calculations  performed  utilized  the  ASA,  we  believe  that  the  trends  will  persist 
in  a  full-potential  calculation,  since  the  volume  of  space  containing  intermediate  spin  po¬ 
larizations  will  be  roughly  the  same.  We  conclude  that  a  good  description  of  the  magnetic 
properties  of  bcc  Fe  is  obtainable  within  the  LSDA,  provided  the  proper  interpolation  is 
used  for  both  the  correlation  energy  density  and  potential. 
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Table  1 


Parameters  for  VWN  interpolation  formula 


parameter 

value 

Aa 

-1/3  n2 

Cct 

[ln(167r/o)  -  3  +  0.531504]/3?r2 

b 

1.13107 

c 

13.0045 

10 

-0.00475840 

a 

(£)* 

12 


Table  2 

Parameters  for  CA  LSDA 


parameter 

Paramagnetic 

Ferromagnetic 

7 

-0.1423 

-0.0843 

k 

1.0529 

1.3981 

02 

0.3334 

0.2611 

A 

0.0311 

0.01555 

B 

-0.048 

-0.0269 

C 

0.0020 

D 

-0.0116 
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Table  3 

Summary  of  LMTO  local  spin  density  calculations  for  bcc  Fe. 


LDA  spin 

interpolation 

lattice  constant 
(a.u.) 

Magnetic  Moment 

(vb) 

Bulk  Modulus 
(Mbar) 

ground 

state 

CA  VWN 

5.3111 

1.98 

2.45 

bcc 

CA  VBH 

5.3706 

1.85 

2.53 

fee 

VBH  VBH 

5.2941 

2.12 

2.52 

fee 

experiment 

5.4169“ 

2.12 

1.68-1.73 

bcc 

a  Reference  22. 
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Table  4 

Summary  of  minimum  total  energy  for  the  ferromagnetic  bcc  and  paramagnetic  fee  calculations. 


LDA 

spin 

interpolation 

Crystal 

Structure 

Energy 

(By) 

CA 

fee 

-44.2974 

CA 

VWN 

bcc 

-44.2993 

CA 

VBH 

bcc 

-44.2945 

VBH 

VBH 

fee 

-44.5952 

VBH 

VBH 

bcc 

-44.5938 

15 


Figure  1 
Figure  2 


Figure  Captions 

.  Total  energy  (minus  the  frozen  core  term)  versus  Wigner-Seitz  radius  for  paramagnetic 
fee  Fe:  (a)  CA  LDA,  and  (b)  VBH  LDA. 

.  Total  energy  (minus  the  frozen  core  term)  versus  Wigner-Seitz  radius  for  ferromagnetic 
bee  Fe:  (a)  CA  LDSA  with  VWN  spin  interpolation,  (b)  CA  LDSA  with  VBH  spin 
interpolation,  and  (c)  VBH  LDSA  with  VBH  spin  interpolation. 
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Energy  (Ry) 


Figure  la 


Figure  1b 
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Figure  2a 


Figure  2b 
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Figure  2c 


